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Effect 
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The dynamical emergence (and subsequent intermittent breakdown) of collective behavior in 
complex systems is described as a non-Poisson renewal process, characterized by a waiting-time 
distribution density ij){r) for the time intervals between successively recorded breakdowns. In the 
intermittent case tp(t) ~ t~^ , with complexity index fi. We show that two systems can exchange 
information through complexity matching and present theoretical and numerical calculations de- 
scribing a system with complexity index fis perturbed by a signal with complexity index fip. The 
analysis focuses on the non-ergodic (non-stationary) case /i < 2 showing that for fis > /ip, the sys- 
tem S statistically inherits the correlation function of the perturbation P. The condition /ip = fig 
is a resonant maximum for correlation information exchange. 

PACS numbers: 05.40.Fb, 05.60.Cd, 02.50.Ey 



The response of complex systems to external excita- 
tions has been one of the central issues of statistical 
physics for the past fiftyyears. These analyses started 
with the work of Kubo [lj and adopted the assumption 
that the unperturbed system can be described by a cor- 
relation function with a finite correlation time. Stochas- 
tic resonance (SR) is one phenomenon that, although 
consistent with the assumptions of statistical physics, 
yields counter-intuitive results Q and has been a popu- 
lar topic of investigation for the past quarter century. SR 
is a strategy for the transmission of information through 
random (Poisson statistics) media that can enhance the 
signal-to- noise ratio Another popular topic of in- 
vestigation has been synchronization of the dynamic el- 
ements in a complex network, whether the elements are 
modeled as linear Q or nonlinear oscillators Q , limit cy- 
cles, or chaotic attractors @. Both the mechanisms of 
synchronization and SR have been used to provide in- 
sight into the behavior of complex networks as, for ex- 
ample, into the dynamics of collections of neurons 
Even more recently the question of how to transmit in- 
formation through a complex network in which the sim- 
plifying assumptions of traditional statistical physics are 
no longer valid has been addressed [1, Q . Barbi et at 
conclude that a complex network described by intermit- 
tent fluctuations with non-Poisson statistics does not re- 
spond to external periodic perturbations [§] in the non- 
stationary (or non-ergodic) regime^ Other investigators 
have reached the same conclusion (Tpj | 

A number of complex phenomena have been shown to 
have non-Poisson statistical properties, including collec- 
tions of blinking quantum dots ll| and collections of 
neurons in the human brain [l2T | . The non-Poisson char- 
acter of the distribution of the sojourn times in the "on" 



and "off' states in BQDs is well known however the 
renewal character of the statistics has only recently been 
established [l3j |. Moreover, Bianco et al. [3| have shown 
using a network of coupled two-state stochastic clocks 
that with the onset of phase synchronization at a criti- 
cal value of the coupling coefficient the dynamics of the 
network becomes that of a non-Poisson renewal process 
operating in the non-ergodic regime. The breakdowns 
of the collective phase structures reset memory of previ- 
ous state changes, consequently yielding a non-Poisson 
renewal process (NRP). NPRs are characterized by a 
waiting-time distribution density between events (colli- 
sions), indicated with ip(r). The regression to equilib- 
rium of an ensemble of NPRs is given by the survival 
probability \t (t) = fl ^){t)cLt. In this letter we focus our 
attention on the case ip(r) ~ with 1 < /i < 2, yield- 
ing inverse power-law relaxation of the form ^(t) w t~ u 
with index v < 1, thereby violating the finite time-scale 
assumption, so often made. One way a non-exponential 
decay of this kind can be expressed is as the sum of in- 
finitely many Poisson components, but if these compo- 
nents are independent, as are the single units in the ab- 
sence of cooperation, no NPR events are generated [lij ]. 
The production of NPR events is a sign of close cooper- 
ation among distinct units, thereby offering a rational as 
to why NPR pro cesses do not respond to periodic per- 
turbations [9l. lioj. as would single independent Poisson 
units. NPR processes reflect a condition shared by the 
phenomenological models of glassy dynamics |l6[, laser 
cooling 
lattices 



17 1 and models of atomic transport in optical 
18[. Other NPR processes have been found at 
the core of the correlations in DNA sequences , heart- 
beat variability [2(| and earthquakes 21 1. 

Herein we extend the investigation of Refs. [{J [l(| to 
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determine how one complex network responds to a per- 
turbation by a second complex network as a function of 
the matching of the measures of complexity of the two 
networks. Here the complex system is an NPR network 
and the measure of complexity is the inverse power-law 
index. More precisely, we consider a NPR system, with 
power-law index /i < 2, and study, as done in the the- 
ory of SR [H, the case where the rate of production of 
jumps, a kind of renewal event, is modulated by an exter- 
nal perturbation. However, we do not use the harmonic 
perturbation of 0, [l(| , but instead use a random signal 
as perturbation. This is to some extent reminiscent of 
the condition of aperiodic SR [22j| . the significant differ- 
ences being that the network satisfies the NPR condition 
with /is < 2 and the perturbation is another NPR pro- 
cess with power index fip < 2. We prove that the NPR 
network is sensitive to the influence of the NPR perturba- 
tion and that the response intensity becomes maximum 
at the matching condition /is = Mp- We refer to this as 
the Complexity Matching Effect (CME), insofar as within 
our theoretical framework, genuinely complex processes 
are NPR systems, fitting the non-stationary condition 
fi < 2: the parameter fj, is the complexity index and two 
systems with the same /j, share the same degree of com- 
plexity. 

Although NPR processes are widespread in the com- 
plexity field, we adopt the BQD picture [IH to illustrate 
the important concepts of jump and collision. The inter- 
mittent BQD fluorescence is a sequence of jumps from 
the "light on" to the "light off" state and vice-versa. We 
use this sequence to generate the dichotomous variable 
£s(i), with the values +1 and —1 corresponding to the 
system in the "light on" (or |+}) and "light off" (or — )) 
state, respectively. These jumps are events whose re- 
newal nature has been proved with compelling arguments 
by the statistical analysis in Refs. [13J. Furthermore, 
the time interval between consecutive jumps is described 
by a histogram with the form of an inverse power law 
^p)(t) ~ with power index fj, < 2. Thus the fluc- 
tuation £s(t) is a NPR process, which, according to Bel 
and Barkai [231 ], violates the ergodic condition. For sim- 
plicity's sake we assume that the distribution densities of 
sojourn times in the two states are identical. The relax- 
ation process from an initial out of equilibrium condition 
is described by 
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T 



t + T 



(1) 



therefore with v = ji — 1. 

This relaxation process is not directly related to the 
distribution density of time intervals between two consec- 
utive jumps, denoted by ip( exp "> (r), but to another distri- 
bution density, indicated by the symbol ip( T )i namely the 
waiting-time distribution for collisions. A collision is a 
non-Poisson event as renewal as a jump, associated with 
a coin tossing that selects the sign of the state produced 



by the collision event. In the absence of perturbation, 
the probability of selecting the state |+) is identical to 
the probability of selecting the state |— ). Consequently, 
the time interval between two consecutive collisions, with 
distribution density ip{r), can be shorter than the time 
interval between two consecutive jumps, but the two dis- 
tributions, i(j(t) and ?/;( ea; P) (r) , have the same power-law 
index \x < 2 [24[. The choice of the form of Eq. ^ 
generates for i/'( r ) the analytical form 



V-(r) = (M - 1) 



{t + TY' 



(2) 



Note that the relaxation function of Eq. (jTJ) , playing 
the role of correlation function, is the survival probabil- 
ity associated with the waiting time distribution of Eq. 
([2|) . Note also that in general the correlation function of 
a NPR is not stationary and "ages" with time, namely, 
with the time lag between preparation and observation. 
In fact, the probability density of the first waiting time 
t\ depends on the time lag between preparation and ob- 
servation. When the observation begins at the same time 
as preparation, the system is said to be "young" and the 
p.d.f. of T\ coincides with ^(r) [24j. The correspond- 
ing survival probability, Eq. |T]), is the young correlation 
function. The distribution density of Eq. ^ is normal- 
ized, and the parameter T, keeping the distribution den- 
sity finite also for r — > 0, gives information on the lapse of 
time necessary to reach the time asymptotic condition at 
which iPs{t) becomes identical to an inverse power law. 
The theoretical results of this Letter are asymptotic in 
time, and do not depend on the specific form adopted for 
V>(t). The two parameters T and \i should be interpreted 
as some short- and long- time properties, respectively. 

In summary, the system S is an ensemble of dichoto- 
mous signals £s(i). Each system of this ensemble is de- 
termined by a sequence of waiting times {t{\ fitting the 
renewal property Pr(j~i,Tj) — Pr(ji)Pr{Tj). This se- 
quence of waiting times generates a series of collisions 
occurring at times {t n } given by t n = J27=o T i- Each 
fluctuation £s (t) is determined by the following prescrip- 
tion: At times t = t n the variable £s(i) gets either the 
value +1 or —1 with equal probability 1/2 and keeps it 
till the next collision. 

To demonstrate CME we assign the role of perturba- 
tion to another NPR of the same kind. Hereafter we 
assign the subscript S to indicate the system properties. 
The system's waiting-time distribution is therefore ips {t) 
of the form ([2]), with parameters T$ and /is. The sub- 
script P is analogously adopted for the perturbation. The 
perturbation is a dichotomous fluctuations £p{t), with 
the values +1 and —1, with the statistics of collisions 
and jumps defined in the same way as for the system 
fluctuations £s(t), thereby involving the survival prob- 
ability \&p(t) and the waiting time distribution ipp(t). 
This perturbation acts on an ensemble of systems char- 
acterized by the survival probability '5 s(t) and the wait- 
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ing time distribution ipg (t) . This perturbation modulates 
the production rate of system jumps with the following 
prescription. When the coupling between system and 
perturbation is switched on, the occurrence of a system 
collision at time t generates the state |+) with the prob- 
ability [1 + e£p(i)]/2 and the state |— ) with probability 
[1 — e£p(t)]/2. The parameter e < 1 is a non-negative 
number defining the interaction strength. If /is = oo, 
namely, when the systems fits the Poisson condition, and 
£p = cos(<jji), this prescription generates the well known 
phenomenon of stochastic resonance 0] . 

Due to the lack of a stationary correlation function, the 
response cannot be ex pres sed by means of the ordinary 
Green-Kubo theory [H.I25I] . Details on how to derive the 
linear response in this case can be found in Ref. [26j |. 
which leads to a theoretical prescription coincident with 
the work of Refs. Thus here we limit ourselves to 

adopting the resulting exact expression that reads 



(&(t)> 



dt'P(t')*s(t-t')Z P (t') 



(3) 



Note that, as in the case of ordinary stochastic resonance, 
we are making an average on infinitely many responses 
to the same perturbation £p(i). Here Pit) denotes the 
rate of collision per unit time, a property of the unper- 
turbed system that remains unchanged even when the 
coupling between system and perturbation is switched 
on. Using renewal theory, it is straightforward to show 
that P(t) — ^2 n=1 i'n(t), w here if> n (t) is the probability 
density of undergoing the n-th collision between times t 
and t + dt 1 which is in turn related to ips(t) through the 
iteration relation i/) n (t) = L ip n -i(t')il)s(t — t')dt' . As a 
result, in the Laplace representation, 



P(u) 



i-i>(u) 



(4) 



In the case /is < 2, the inverse Laplace transform of Eq. 
(g]) yields, for t — > oo, P(t) oc l/£ 2 ~ Ms , indicating that 
the rate of collision production decreases as a function 
of time, which is a different way of expressing the aging 
condition mentioned earlier. 

CME implies that the perturbation is also an erratic 
signal. Thus, the only possible way to make an analyti- 
cal prediction rests on using infinitely many realizations 
of the same perturbation signal £p(t) so as to produce 
a Gibbs ensemble of perturbations P as well as a Gibbs 
ensemble of systems S. Let us assume that all the pertur- 
bation signals £p(i) are properly prepared at time t = 
in the condition £p(0) = 1, while the system is prepared 
in a condition with half the systems in state |+) and half 
in the state |— ), so that (£s(0)) = 0. In practice, as 
an effect of this prescription, we average Eq. over 
all possible realizations of the perturbation £p (t) , which 
yields the double average ((£s(i))) in terms of the simple 



average (£p(i)). It is known [24[ that (£p(i)) = *p(t). 
Thus, due to the linear nature of Eq. ^ we obtain 



«&(*)» 



dt'P(t')^ s (t-t')^p(t') 



(5) 



We are now properly equipped to make analytical pre- 
dictions. Due to the time convolution nature of Eq. ([5]) 
it is convenient to evaluate its Laplace transform, and to 
consider the asymptotic limit u — ► 0, which by inverse 
Laplace transform yields information in the time asymp- 
totic limit t — v oo. The calculations are done using the 
approach illustrated in Ref. [27j and yield, for t — > oo, 



«&(*)» = 



(6) 



We see that the two terms of the r. h. s. of © are re- 
spectively proportional to the asymptotic tails of \&5 and 
^ p with amplitudes given by the analytical expressions: 



As 
A P 



Ms) 



r( MP )r(i- Ms ) 
r(i — /xp + fj, s ) 



and 



r(^s - tip) 



r(2 - Mp)r(Ms) r(/i S )r(i- M + p) 



(7) 
(8) 



As and Ap have opposite signs, with the former being 
positive for jis < fJ-p and negative for fj,g > fip. In both 
cases the dominating term is the one with the positive 
pre-factor, signaling that the system is able to inherit the 
statistical properties of the perturbing system when //g > 
/ip, while maintaining its own statistical properties when 
Ms < Mp- At /ig = /ip both terms diverge and this is due, 
as we shall see by means of a numerical calculation, to 
the critical slowing down corresponding to the matching 
condition. 

To shed light on the physical meaning of this analytical 
result, it is convenient to compare Eq.® to the exact 
solution of Eq. ([S]) in the case when the system satisfies 
the Poisson condition with rate gs and the perturbation 
is another Poisson system with rate gp. In this case 



(9) 



As = gs/(gp-gs), A P = -A s , * s (t) = exp(-g s t) 
and ^>p(t) = exp(-gpt). We see that the excitation- 
relaxation process corresponds to the superposition of 
two exponential contributions with coefficients of equal 
amplitude and opposite sign, with a diverging maximum 
at g P = g s . This example, closely connected to the ape- 
riodic stochastic resonance [22| , shows that the maximum 
response of the system S to the excitation £p(i) is ob- 
tained when the rate of the perturbation is identical to 
the rate of the system. This suggests the adoption of 
Eq. ((9]) as a fitting function of the Monte Carlo calcu- 
lation that we use to simulate CME. In this latter case 
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^s{t) and ^p(t) are given by Eq. dXJ) , with subscripts 
S and P, respectively. Using Eq. ©, we fit the values 
of As and to the Monte Carlo data. CME is ade- 
quately described by Fig. 1, which is obtained by keeping 
Us = 1-6 and moving /jp from [Lp = 1.1 to [is — 1.95. 
We plot the intensity of the slowest response contribution 
as a function of \xp. We see that there is an impressive 
similarity with the Poisson resonance phenomenon. The 
maximum response is obtained by realizing the matching 
condition /ip = [as, in analogy with the resonance con- 
dition gp — gs- Two fitting examples, [is = 1-35 and 
Us = 1.85 are shown in the inset of Fig. 1. 
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FIG. 1: Inset: fitting of Eq. Q (solid lines) to Monte Carlo 
data (open circles) using Ts = Tp = 1, /j,s = 1.6 with fip = 
1.35 (upper) and [tp = 1.85 (lower). Dashed lines are the 
asymptotic dominant term in Eq. <(9j . Our Monte Carlo used 
10 7 system-perturbation pairs. Main figure: Amplitudes Ap 
(squares), As (triangles), Eqs. JH} (solid line) and (J7| (dashed 
line) as a function of [ip, with [is = 1.6. 

The CME is a general property of a large class of com- 
plex systems and is not limited to the intermittent flu- 
orescence used to motivate our discussion. For exam- 
ple, other NPRs are given by return times for a random 
walk, either in regular lattices or in a complex network, 
regardless of whether the diffusion process is normal or 
Fractional Brownian Motion 28]. In general, events are 
generated by a slow motion, which reflects the creation 
and the subsequent dynamics of self-organized structures. 
Events either break these collective modes or reflect the 
system entering a regime with high rates of entropy in- 
crease. This can be, for instance, a passage in a chaotic 
sea, the effect of a coin-tossing or a catastrophic " quake" 
in the system structure Whatever the case, the sys- 
tem memory is reset, while the non- Poisson behavior of 
■0(t), with uncorrelated waiting times, is responsible for 
long-range correlations in the observable sequence. 
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